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ABSTRACT 

We present new X-ray and Far Ultraviolet (FUV) observations of T Tauri stars 
covering the age range 1 to 10 Myr. Our goals are to observationally constrain the 
intensity of radiation fields responsible for evaporating gas from the circumstellar disk 
and to assess the feasibility of current photoevaporation models, focusing on X-ray 
and UV radiation. We greatly increase the number of 7-10 Myr old T Tauri stars 
observed in the X-rays by including observations of the well populated 25 Ori aggregate 
in the Orion OBla subassociation. With these new 7-10 Myr objects, we confirm that 
X-ray emission remains constant from 1-10 Myr. We also show, for the first time, 
observational evidence for the evolution of FUV radiation fields with a sample of 56 
accreting and non-accreting young stars spanning 1 Myr to 1 Gyr. We find that the 
FUV emission decreases on timescales consistent with the decline of accretion in classical 
T Tauri stars until reaching the chromospheric level in weak T Tauri stars and debris 
disks. Overall, we find that the observed strength of high energy radiation is consistent 
with that required by photoevaporation models to dissipate the disks in timescales of 
approximately 10 Myr. Finally, we find that the high energy fields that affect gas 
evolution are not similarly affecting dust evolution; in particular, we find that disks 
with inner clearings, the transitional disks, have similar levels of FUV emission as full 
disks. 
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Introduction 



Circumstellar disks, formed as a consequence of the star formation process, evolve from an 
initially gas rich phase to a system composed of small, planetary sized bodies, with little remaining 
gas. The physical processes responsible for depleting circumstellar gas from the initial minimum 
mass solar nebula (MMSN) to the low quantities remaining in debris disks, th ose with significant 



dust evolution and gas depletion, have been the focus of gas dispersal models (|Clarke et alJ 12001 



Alexander et al. 



2006 



Gorti &: Hollenbach 



2009a 



Gorti et al 



2009b 



Owen et al 



Models for 



the evolution of the dust (jPullemond k Dominikll2005l : lD'Alessio et al.ll2006l : lRobitaille et alJl2007l ) 
have been successfu l ly utilized in the interpretation of a variety of Infrared ( I R) spectral shapes 
dSargent et alJbood : iHernandez et alJboOTal : lEspaillat etaDbooi ICurrielbood : IPurlan et allbood : 
Najita et al.ll20ld ). These models indicate that while some of the material is locked into planets 
and still more is accreted onto the star through viscous evolution, it is likely that photoevaporation 
is the most significant process playing a role in the dispersal of the gas disk. There are several 
contributors to photoevaporation, including X-rays or Extreme Ultraviolet (EUV) photons from 
the central star and Far Ultraviolet (FUV) e mission produced i n an accretion sho ck at the stellar 
surface of classical T Tauri stars, or CTTS dOwen et al.ll2010l : iGorti et al.ll2009bl ). At any given 
time, contributions from several different photoevaporation mechanisms will occur, but the overall 
timescales when each contribution is greatest may vary. For example, accretion related radiation 
fields will decay as the mass accretion rate (M) decreases. Alternatively, X-ray emission may 
increase as the accretion related radiation fields drop due to increased magnetic activity in non- 
accreting T Tauri stars (WTTS), observed as a possible increase in X-ray luminosity from 1-100 



Myr in the Chandra Orion Ultradeep Project (COUP) survey of young stars (jPreibisch &: Feigelson 



20051 ). The changing environment may play a role, as FUV and X-rays dominate photoevaporation 
early at M < 10~^ Mq yr~^, while EUV photons are absorbed in the accretion streams and do not 
become a significant contributor until M drops to < 10~^ M0yr~^ ( Gorti IIollenbachll2009al ). 



With the large number of factors to consider in gas disk evolution, the intensity of UV and X-ray 
radiation fields and the timescales relevant to their evolution are essential parameters in the theory 
of photoevaporation. Not only observations of the radiation fields, but observations of the frequency 
of circumstellar disks place constraints on disk evolution theory, with only 10% of sources retaining 
disks at 10 Myr (IHernandez et al.ll2007bl). Full (or primordial) disks and debris disks are observed 
at both 5 and 10 Myr (jHernandez et al.ll2006l . l2009l ). indicating that disk dispersal mechanisms must 
occur on a number of timescales in a single star forming region. Transitional and pre-transitional 
disks, which due to the presence of a gap are interpreted as an interme diary stage between full and 



debris disks (IStrom et al.lll989l : ICalvet et al.ll2002l : 



Espaillat et al.ll2008l ) , are observed to have a low 



frequency indicating t hat when disk dissipat i on begins it occurs qu ickly, taking as little as 0.5 Myr 
(jLuhman et al.l l2010l : IHernandez et al.l l20ld : iMuzerolle et al.l l2010l ) . Alternatively, disk evolution 
may proceed along another path in which a similar decline in IR fluxes from all disk radii iridicate s 
that no gap or hole is formed (jSicilia-Aguilar et al.l l2010l : iManojl l2010l : IHernandez et al.l |2010| ) . 
With the diversity observed in disk evolution, different dispersal mechanisms may be necessary. 
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In this paper we investigate the relevant timescales for the evolution of disk dispersing emission 
from T Tauri stars with new X-ray and FUV observations spanning 1-10 Myr, and include previous 
observations of a larger sample of observations extending to young stars at 1 Gyr. 

While EUV emission is almost impossible to observe due to interstellar hydrogen absorption, X- 
ray observations are numerous. X-ray emission in T Tauri stars is significantly stronger, up to 1000 
times, than the X-ray emission from main seq uence stars, and come s from hot plasma with temper- 
atures up to 60 MK or energies up to 5 keV (jFeigelson et al.ll2002l ). The current paradigm is that 
this hot plasma is produced through enhanced magnetic activity, wit h stronger magiietic re connec- 
tion events the likely cause of observed short period flaring activity (jFeigelson et al.lll993l ). There 
is some evidence that strong soft X-ray emission is produced in an accretion shock on the stellar 
surface (iGiinther et al.ll2007l ). Currently observations of a soft X-ray component which dominates 
over the coronal emission are unique to the star TW Hya, although ot her sources have exhibite d 
accretion produced X-ray emission, though at a significantly lower level (IRobrade &: Schmittll2006l ). 
Surveys of young, <10 Myr star forming regions have shown no evidence of a decline in the X-ray 
emission as M decreases, indicating that accretion relat ed X-ray emission is not a s i gnificant con- 



tribu tor to the total X-ray flux in most circumstances (jPreibisch &: Feigelson 



2005 



Mercer et al 



20091 ). While roughly constant between 1 and 10 Myr, eventually the X-ray emission does decrease. 
When this decline begins has not been well characterized so far because the only 10 Myr o ld X-ray 
sources studied were the eleven sources in the TW Hydra Association (jWebb et al.lll999l . TWA). 
Previously, the cluster nearest i n age to TWA with well c haracterized X-ray properties was the 
6-8 Myr -q Chamaeleon cluster (jLopez-Santiago et al.ll2010l ). with fourteen sources. In this paper 
we analyze X-r ay observations frorn another older cluster, the 25 Ori aggregate in the Orion OBla 
subassociation (iBriceho et al.l 120071 ). At 7-10 Myr, 25 Ori still has 4-5 remaining CTTS and a 
large number of WTTS, making it the most populous ~10 Myr old star forming region known. 
Extinction in this ofi^-cloud area is low, making it an ideal region for studying X-rays from T Tauri 
stars in older, evolved populations. 

Less is known about FUV emission from T Tauri stars than X-ray emission because, until 
recent instruments like the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope 
(HST), most FUV detectors lacked the sensitivity to observe dim sources. Large FUV samples have 
been obtained with the Inte rnational Ultraviolet Explorer (lUE) throug h which FUV line emission 
has been well characterized (IValenti et al.ll2000l : iJohns-Krull et al.ll2000l ). but little to no signal was 
achieved in the continuum at A < 1750 A for the WTTS. Space Telescope Imaging Spectrograph 
(STIS) observations of T Tauri stars have also been analyzed; however, long exposure times were 
required for all but the brightest sources, lira iting the ability of STIS to observe a large population 
(jCalvet et ahlbood : Eerczeg et ahlbood . bood ) . Even though lUE and STIS observations are hmited, 
the T Tauri properties observed in these high resolution spectra were used to interpret low resolution 
spectra obtained in a large sample of CTTS, WTTS, and debris disks using the Advanced Camera 
for Surveys Solar Blind Channel (ACS/SBC) prism on HST. While low in resolution, with its 
excellent sensitivity ACS/SBC obtained some of the first observations of WTTS in the FUV with 
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the continuum detected from 1230 - 1800 A (jlnglebv et alibood ). 



The origin of FUV emission in T Tauri stars has contributions from s everal sources. In Ingleby 
et al. (2011), we show that in addition to an accretion shock component (jCalvet et al.ll2004l'). Hq in 



the inner region of CTTS disks contributes to the FUV Une and co i itinuum emission (|Abgrall et al 



19971 : llnglebv et al.l l2009l : iBergin et alJ l2004l : iHerczeg et al.l 12004 I2OO2I : iFrance et al 



201 



I. 



The 

evolution of FUV emission has not been well studied, as FUV spectra of T Tauri stars greater than 
1-3 Myr are rare, and previously included only four T Tauri stars in the TWA region and older 
debris disks (jlnglebv et al.ll2009l : iHerczeg et al.ll2004l ). 



In this paper, we discuss the evolution of high energy radiation fields by analyzing the X-ray 
and FUV luminosities of a large sample of sources (the largest FUV sample to date) with ages 
between 1 Myr and 10 Myr, including a small sample of post T Tauri objects with ages up to 1 
Gyr. In ^we discuss the new X-ray, FUV, and optical observations used in our analysis. In ^3.11 
we present new X-ray selected T Tauri stars in 25 Ori and Chamaeleon I. In ^3.2l we present results 
of spectral analysis performed on high count rate X-ray sources and we discuss the calculation of 
X-ray and FUV luminosities in ^3.31 and ^3.4[ Finally, in iQlwe analyze the high energy radiation 
fields and discuss whether the observed evolution is consistent with current disk evolution theories. 



2. Source Sample and Observations 
2.1. X-ray Data 



Chandra observations of four fields in the 7-10 Myr 25 Ori star forming association were ob- 
tained in 2008 January with ACIS-I and S chips. Observations were done in FAINT mode with 10 
ks exposures in GO Program 09200907 (PI: Hartmann). Additionally, observations of two fields in 
the 2.5 Myr Chamaeleon I star forming region were obtained in 2009 February with ACIS-I and S 
chips. Again, the observations were taken in FAINT mode in GO Program 09200909 (PI: Calvet) 
with 25 ks exposure times. We analyzed the evt2 files provided from the Chandra version 7.6.11.4 
processing pipeline. X -ray detections were identified using the Chandra Interactive Analysis of 
Observations software ([Fruscione et al.l 120061 . CIAO) script WAVDETECT with t he defaul t signi f- 
icance thresho l d set to 10~^. We cross-correlated the X-ray dete ctions with the iLuhmanI (120041 ). 
Luhman et"al] hooi ). and the CIDA VariabiUty Survey of Orion (jBriceho et ahlboosl . CVSO) op- 
tical catalogs and identified previously studied T Tauri sources. 

In Chamaeleon I, the selected fiel ds focused on tw o CTTS with disks showing evidence for 
inner clearing, SZ Cha (T6) and T35 (jKim et al.ll2009l ). We detected 15 known T Tauri stars in 
the two fields (out of ~60 X-ray sources identified by WAVDETECT), five of which are accretors 
(jLuhmanll2004l ) . The sources are listed in Tableland their positions in the H-R diagram are shown 
in Figure m based on spectral typ es and lumi n ositie s from the literature (see references in Tabled]). 
Using the evolutionary tracks of ISiess et al.l ()200d ). we determine that the sample covers a large 
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range of T Tauri masses from a 1-1.5 Mq upper limit down to ali nost substellar masses ^ with a 
lower limit of 0.1-0.2 Mq. This mass range is similar to that used in lPreibisch &: FeigelsonI (|2005l ). 
in which the evolution of X-ray emission in the low mass population of the COUP sources was 
analyzed. 

Each of the four fields in th e 25 Ori aggregate was centered on a know n CTTS, specifically 
CVSO 35, CVSO 206, OBla 1630 dBriceho et ahlbood : iHernandez et al.ll2007al l. and the transitional 
disk CVSO 224 (jEspaillat et al.1 [200^). Of ~150 X-ray sources located by WAVDETECT, we 
id entified 36 prev i ously known T Tauri stars, in addition to the 4 targeted sou r ces, fro m the catalogs 
of briceho et al.1 (|2005l . boO?! ). Briceho et al. (2010), and iHernandez et al.l (|2007al ) for a total of 
40 7-10 Myr classical and weak T Tauri stars. Of the remaining unidentified X-ray detections, 
six are possible newly detected T Tauri stars (see §2.31 for further discussion). The locations of 
the 25 Ori sources in the H-R diagram are shown in Figure [2j Sources pub lished in 



(120051) had know i i spect ral types and luminosities but sources published in iBriceho et al 



Calvet et al 
2007l y 



Hernandez et al.l (j2007al ). and Briceho et al. (2010) had only spectral types and optical colors. To 



get luminosities for thos e sources, we first used V — I colors to get the reddening (Ay) assuming 



Survey (jSkrutskie et al 



photosp heric colors froinlKenvon &: HartmannI (|l995l ) . With Ay , we corrected Two Micron All Sky 



20061. 2MASS) J magnitudes for reddening usi ng the Mathis extinction law 



with Rv = 3.1 (jMathisI Il990l ) and calculated the bolometric fiuxes (jKenvon fc HartmannI Il995l ). 
Our sample of 25 Ori sources is similar to the CTTS masses of Taurus, from ~0. 2-1.0 Mq. Tabled] 
lists the complete set of sources detected in the Chamaeleon I and 25 Ori fields which are discussed 
in this paper, along with published spectral types and disk classifications. 

After matching X-ray detections to previously identified T Tauri stars, we used the CI AO script 
DMEXTRACT to obtain the counts within the extraction regions determined by WAVDETECT 
for each known T Tauri source. We included a background region defined as a 5" annulus around 
each extraction region, allowing us to determine the net counts and calculate the poisson errors 
which we used to estimate uncertainties in the X-ray luminosities ( ^3.3p . Of the 15 Chamaeleon 
sources, 5 had greater than 250 counts and the spectra of these sources were extracted using the 
CIAO script PSEXTRACT, using the 3a detection regions determined in WAVDETECT, and the 
5" background annuli. Each extracted spectrum was re-binned to 20 counts per bin using the 
GRPPHA task in FTOOLS. The re-binned spectra were then analyzed in the X-ray spectral fitting 
package, XSPEC (jArnaudlll996l ). Only three of the 25 Ori sources had counts in excess of 250, 
likely due to the shorter exposure times and greater distance to the sources, and were analyzed 
with PSEXTRACT and XSPEC. 



2.2. UV Data 

We present new FUV observations of ten sources in the 2.5 Myr Chamaeleon I and two in the 
4 Myr Chamaeleon II star forming regions obtained with the PR130L grating on the Advanced 
Camera for Surveys (ACS) Solar Blind Channel (SBC). These observations were completed be- 
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tween 2008 December and 2009 February in GO Program 11145 (PI: Calvet). The FUV coverage 
extends from 1230-1800 A with R~300 at 1230 A and ~50 at 1700 A. The Chandra and ACS/SBC 
Chamaeleon I samples have some overlap but there are sources unique to each sample. We show 
ah the Chamaeleon I and II sources in Figure [H identifying the X-ray and FUV sources and those 
included in both samples. With ACS/SBC we observed sources from 0.3-2.0 Mq in Chamaeleon I 
and approximately solar mass sources in Chamaeleon II. 

We combine these new observations with pr e viously observed AC S/SB C and STIS s a mples 
of T Tauri stars analyzed in llnglebv et al.1 pOOOl ). ICalvet et al.l (j200J), and iHerczeg et al.l (|2002l . 
2004 ) . The total sample of sources with FUV spectra consists of 56 sources, eleven of which are 
WTTS or debris disks, and these are listed in Tabled! with sources observed by STIS distinguished 
in bold font. 



2.3. Optical Data 



A large number of X-ray detections corresp onded to unidentified counterparts in IR and optical 
catalogs. Using 2MASS and CVSO photometry (IBricerio et al.ll2005l ). we produced color-magnitude 
diagrams in V and V — J to select sources which fell above the main sequence. In Chamaeleon I, 
one such object met the color requirements (K. Luhman, personal communication) and in 25 Ori 
there were six possible new T Tauri stars detected. 

For three of the six T Tauri candidates in 25 Ori, we followed up the X-ray observations by 
obtaining optical spectra wit h the Magellan Echel lete (MagE) on the Magellan- Clay telescope at 
Las Campanas Observatory (IMarshall et al.l l2008l ) . MagE provides R~4100 spectra from 3200- 
9800 A, including the Li 6707 line, for confirmation that the sources are young, and Ha necessary 
to classify them as T Tauri stars. In addition to the MagE spectra of new T Tauri stars, we 
obtained R~ 35,000 MIKE (The Magellan Inamori Kyocera Echelle) observations of 16 sources in 
our sample, primarily of accretors in the 25 Ori and Chamaeleon I and II regions. 



The optical data were reduced using the Image Reduction and Analysis Facility (jTodylll993l . 
IRAF) tasks "ccdproc", "apflatten" and "doecslit". Bias corrections were completed with "ccd- 
proc" . Flats were created using "apflatten" , in which the orders were traced using observations of a 
quartz lamp for the MIKE flat frames and the 10 reddest orders of MagE. For the five bluest orders 
of MagE, we observed a Xe lamp with a wide slit (5" as opposed to the science slit at 1") and out of 
focus for the trace. The task "doecslit" was used to extract the science spectrum, extinction correct 
for airmass, subtract the background, and dispersion correct the spectra using a ThAr comparison 
lamp. The blaze function was fit using a bright O or B standard star and used to correct for the 
blaze in the science exposures. 
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3. Analysis of Observations 
3.1. Results from Optical Observations 

MIKE optical spectra were used to confirm status of WTTS or CT TS and obta i n acc retion 



rates for the CTTS using the relation between M and Ha 10% width of lNatta et alj (|2004l ). We 
show Ha line profiles in Figure [3] for the accreting sample of MIKE sources. The relation between 
Ha width and M is more appropriate for low values of the mass accretion rate; more precise values 
of the accretion rate will be obtained in a forthcoming analysis of the optical veiling of these sources 
(Ingleby et al. 2011). 

The three candidate members of 25 Ori with MagE optical spectra, discussed in ^2.31 showed 
Li absorption and but only two (2MASS J05243490+0154207 and 2MASS J05270173+0139157) 
had weak Ha emission, indicating that they are newly identified WTTS. The third source (2MASS 
J05244265+0154116) has Ha in absorption but may be a young member of 25 Ori still, due to its 
early spectral type; additional analysis is needed for confirmation. This final source is not included 
in any further analysis of the 25 Ori X-ray observations. The Li and Ha equivalent widths, along 
with the derived spectral types are listed in Table [2] for those sources which have optical spectra. 
All of the candidate X-ray selected T Tauri stars, in both Chamaeleon I and 25 Ori, are listed in 
Table [2j MagE spectra will be obtained for the remaining sources in need of optical spectra. 



3.2. X-ray Spectral Analysis 



Eight sources in our Chamaeleon 1 and 25 Ori samples had strong X-ray fluxes allowing for 
spectral analysis to determine temperatures for the emitting plasmas. Three of th ese sources are 
accreting (CR Ch a, SZ Cha, and TW Cha) and the remaining five sources are not (|Luhmanll2004l : 
Calvet et al.ll2005l ). With this sample of high count sources, we looked for any observable distinc- 
tion between the CTTS and WTTS. In XSPEC, each source spectrum was fit using either one or 
two APEC (Astrophysical Plasma Emission Code) models characterized by different temperatures; 
the choice of models was motivated by th e two-temperature plasma observed in X-ray emission 
from Pleiades sources (jBriggs &: Pvd 120031 ). Assuming a model which includes multiple temper- 
atures also allows for the distinctio n of a soft accretion related emission component, if present. 



similar to that observed in TW Hya (jCiinther et al.l 120071 ). We included an absorption component 
using the XSPEC model, PHABS, with the equivalent neutral hydrogen column densitjQ obtained 
from the HEASARC (High Energy Astrophysi cs Science Archive Research Center) Nh calculator 



(Kalberla et al 



2005 



Dickey &: Lockman 



1990|). 



*The equivalent neutral hydrogen column density (defined as Nh) is the column density of hydrogen required to 
give the same opacity as that resulting from the true composition of the absorbing gas, which has contributions from 
He, C, N, O, and Ne in addition to less abundant gasses. 
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The models with the best statistical fit to the observations (or lowest Xred) given in Table 
[31 The best fit models for CR Cha, SY Cha, SZ Cha, and TW Cha produced good Xred^^ but the 
remaining sources had either poor fits or model properties discrepant from the larger sample. For 
CHXR 30A, when the absorption model was for ced to have a column density of Nh = 8 x 10^*^ cm~^ 



(Kalberla et al 



2005 



Dickey Lockmanlll990l ) , a 64 keV plasma was necessary to f it the spectrum 



great er than ten times higher than the highest expected plasma temperatures (jFeigelson et al 



2002)- Allowing the column density to vary, we found a best fit model with Nh = 12.2 x 10^^ cm 



fifteen times higher than Nh in the four Chamaeleon sources with good fits. While high, a greater 
value for Nh is consistent wi th CHXR 30A being extincted by Aj ~ 3 magnitudes while the other 
sources have Aj = 0.0 - 0.5 dLuhmanlboOTl '). The APEC model fits to CVSO 38 res ulted in a large 



temp erature for the hot component, a temperature typical of flares in T Tauri stars (jGetman et al 



20081 ). indicating that the observations may have caught it while flaring. The APEC components 



and best model fit to the highest count rate X-ray source (CR Cha) are shown in Figure [H For 
CR Cha we find that a two-temperature model improves the Xred ^ factor of 3 compared to a 
single temperature model. Like CR Cha, the models of the four sources with the best xied 
each comprised of a cool, 0.3-0.9 keV, component and a hotter 2-5 keV component. We note that 
the hottest emission (excluding the possible flare component from CVSO 38) comes from SZ Cha 
which, as a G2, is the source with the earliest spectral type. 

Previous analyses of X-ray spectra from CTTS have shown that the presence of accretion alters 
the X-ray spectrum by producing emission lines fo rmed in the high densities of the accretion shock 
(jRobrade Schmittll2006l : iBrickhquse et al.ll201C). Also, cool X-ray components with kT ~0.2-0.3 



have been attributed to accretion ( Giinther et al 



20071 : iRobrade fc SchmittJl200fil ). We see no clear 



evidence of accretion related soft X-ray emission in CR Cha, SZ Cha, and TW Cha, the known 
accretors in our sample. The cool APEC components are ~2-10 times hotter than the temper- 
atures expected of plasmas heated in an accretion shock, even when considering material falling 
faster tha n the free fall veloc it y which is necessary to produce the soft X-ray shock emission in 
TW Hya (jCiinther et al.l 120071 ). iBrickhouse et al.l (|2010l ) obtained high resolution Chandra spectra 
of TW Hya and identified line emission which is attributed to accretion. Additional observations 
obtained with XMM- Newton of a sample of four CTTS, including both TW Hya and CR Cha 
(which we analyze in Figured]), showed that shock formed line emission may be needed to confirm 
accretion for sources unlike TW Hya, where the soft accretion component does not strongly exceed 
the coronal emission. The resolution of our ACIS spectra prohibited accretion line analysis, which 
would be informative, particularly for TW Cha, the accretor with a cool component most similar 
to that in TW Hya. Interestingly, the sources with the coolest components (SY Cha, 2MASS 
J05244498+0159465 and RX J0526+0143) are WTTS and have been identified as such by obser- 
vations of Ha, indicating that a detectio n of a cool component alo ne cannot confirm accretion and 
likely line detections are needed as well (jRobrade fc Schmittll2006l ). 



Instead, our results are consist ent with the two-temperature coronal emission observed in 
Pleiades sources (jBriggs &: Pyell2003l ): however, for Chamaeleon, the temperatures of the hot com- 
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ponent are approximately 2 times higher than the hot components in the Pleiades sources. This 
may be indicative of theori zed cooling coronal gas, producing a softer spectrum as X-ray sources 
evolve (jKastner et al.l 120031 ) . The results of our spectral analysis show that we can use the same 
treatment to find X-ray luminosities for both the WTTS and CTTS since we see no evidence for 
an accretion related component. 



3.3. X-ray Luminosities 



The majority of the X-ray sources had <250 counts, so instead of attempting detailed spec- 
tral analysis on dim sources we assume typical X-ray temperatures as determined for the sources 
discussed in ^3.2[ We assume there is no discernible difference between the X-ray properties of the 
CTTS and WTTS in the spectral region observed by ACIS, as indicated by the results shown in 
Table [3l Using WebPIMMS on the HEASARC website we calculated unabsorbed fluxes between 
0.2 and 10 keV for the entire sample, correc ting for Nh = 0.8 x 1 0^^ and Nh = 1.1 x 10^^ cm~^ 
for Chamaeleon I and 25 Ori, respectively (iKalberla et al.l l2005l : iDickey &: LockmanI Il990l ). We 



assumed a single temperature APEC model with kT=2.2 keV and solar abundance. From the 
unabsorbed fluxes. X-ray luminosities were calculated assuming distances of 160 and 330 pc for 
Chamaeleon I and 25 Ori, respectively; the results are listed in Table [TJ For the sources in Table El 
which were fit with XSPEC, we find that the when using WebPIMMS to find fluxes, the discrep- 
ancy is a factor of 2-3. The distribution of unabsorbed X-ray luminosities is shown in Figure \5\ 
which indicates that our Chamaeleon I sample is not complet e. To increase the sample, we include 
X-ray obse rvations of this region from iFeigelson et al.l (|l993l ) in our analysis. The ROSAT count 
rates from IFeigelson et al.l (|l993l ) were converted to ACIS-S count rates using WebPIMMS before 
luminosities were calculated. 

In addition to the X-ray luminosities calculated here, we include X-ray luminosities from the 
literature for sources with FUV spectra which were not observed with Chandra, primarily in Taurus, 
TWA, and nearby debris disks (shown in italics in Table [1]). X-ray fluxes of Taurus sources were 



calcu l ated as part of the XMM Newton Extended Survey of Taurus (XEST) project (jCiidel et al 



2007 


)• 


Giidel et al. 


(2007) 



converted into Chandra ACIS-S count rates using WebPIMMS, assuming an input XMM-Newton 
energy range of 0.4-10 keV and output Chandra energy range of 0.2-10 keV, and using the same 
model parameters as were used to cal culate the fluxes fro m our Chandra observations. X-ray fluxes 



for sources in Taurus not provided by I Giidel et al.l (120071 ). along with sources in Ch amaeleon I and 



II, TWA, and field deb r is disk s were converted from ROSAT PSPC count rates (jFeigelson et al 



1993 



Voges et al.l Il999l . l2000l ) into Chandra ACIS-S count rates, again using WebPIMMS but 
using the softer input energy of ROSAT (0.1-2 keV). Once converted to ACIS-S count rates, the 
unabsorbed X-ray fluxes were calculated as discussed for our new Chandra observations and the 
luminosities are given in Table [TJ We also include X-ray luminosities from ot her young clusters i n 
our analysis of X-ray evolution; those from the ~ 4 Myr cluster, Trumpler 37 (jMercer et al.ll2009l ). 
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and the 6-8 Myr r] Cha region (jLopez-Santiago et al.ll20ld ) 



3.4. FUV Luminosities 



Our FUV sample contains sources in Taurus, the Chamaeleon I and II clouds, 25 Ori, TWA, 
objects in the outskirts of the Orion Nebular Cluster (ONC), and nearby WTTS or debris disks. 
We calculated FUV luminosities for these sources between 1230 and 1800 A after degrading the 
resolution of the STIS spectra to that of ACS / SBC for c omparison. S p ectra were corrected for 
reddening using j4y's provided in llngleby et al.l ([20091) and lCalvet et al.l ([200J) for Taurus, ONC, 
25 Ori, and the nearby WTTS and deb ris disks. For the Cham aeleon I and II sources, Ay's were 
calculated from published photometry ([Gauvin &: Stromlll992l ). For sources in Taurus, the ONC, 
and Chamaeleon I and II we corrected for reddening using the extinction curve observed along the 
line of sight towards HP 29 647 with Ry = 3.6, appropriate for regions with associated molecular 
clouds (jWhittet et al.ll2004l ). For sources in the dispersed 25 Ori association, we used the Mathis 
extinction law with Ry = 3.1 (jMathislll990l ). We note that UV extinction laws, especially in the 
FUV, are very uncertain and introduce the largest error in our FUV luminosity calculations. We 
assume a ±0.2 magnitude error in Ay, which dominates other sources of error. For four Chamaeleon 
sources, we failed to detect the FUV emission and upper limits are provided for these sources. We 
calculated the flux in the detector noise and de-reddened the noise by the calculated Ay's. In 
this manner, we are accounting for any FUV emission washed out by extinction. Calculated FUV 
luminosities are given in Table [H with upper limits on the luminosities provided when appropriate. 



4. Implications of Observations 
4.1. Evolution of X-ray and UV Emission 



The intensity of high energy emission from T Tauri stars has important implications for the 
dispersion of circumstellar disks and the planets which may be forming within them. Over a 
time period of 1-10 Myr, during which these processes are expected to occur, we investigate the 
evolution of the X-ray and FUV radiation fields. Figure [6] shows the median and full range of X-ray 
luminosities for the sources in each star forming region with spectral type K5 or later. The sample 
of Tr37 may include sources with earlier spectral types (early K); however, spectral types are not 
well established for this sample. In addition, it was confirmed that several of the Tr37 sources were 
flaring during the observations (jMercer et al.ll2009l ). contributing to the high X-ray luminosities. 
Figure [6] reveals that there is no decrease with age in either total X-ray luminosity or the X-ray 
luminosity normalized by the stellar luminosity, up to 10 Myr. Most importantly, our addition of the 
7-10 Myr 25 Ori X-ray observations supports previous cl aims that X-ray radiation from the central 



star r emains strong throughout the disk dispersal phase ([Mercer et al 



200S 



Preibisch &: Feigelson 



20051 ). While there does appear to be a slight increase in the normalized X-ray luminosity with 
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age, the trend is statistically insignificant, with a Pearson correlation coefficient of only 0.2. 



Figure [7] shows the evolution of FUV luminosities from 1 Myr T Tauri stars to 1 Gyr young 
stars and debris disks. While there is an observed decrease in FUV luminosity out to 1 Gyr, a 
decrease with age can also be seen among the 1-10 Myr sources, with the WTTS characterized by 
the lowest FUV luminosities. By the WTTS phase, FUV luminosities are expected to be entirely 
chromospheric. With no additional accretion, FUV emission is no longer produced in an accretion 
shock and because the end of accretion indicates the depl etion of inner disk gas, FUV emission from 
hot molecular gas no longer contributes to the FUV flux (jlnglebv et al.ll2009l : iPascucci et al.ll2006l ). 
In addition to the decline with age, a correlation between FUV luminosity and Lace is observed, 
with a Pearson correlation coefficient of 0.78 (calculated with the exclusion of sources with only 
upper limits), indicating that FUV emission is strongly correlated with Lace- Values of Lace were 
obtained from the literature (see references in Table [T]) or calculated from Ha 10% line widths 
(|Natta et ahlbooi ^3T]) . 



This correlation between FUV luminosity and Lace is not surprising, given results that FUV 
emission is produced in t he hot gas of the accretion shock (jCalvet et al.l 12004 : ICalvet &: Gullbring 
19981 : iHerczeg et al.ll2002l ). Therefore, the decrease in FUV lum inosity between 1 and 10 Myr can be 



understood in terms of the dec r ease of M in this time period (jHartmann et al 



2005 



1998 



Gal vet et al 



Siciha-Aguilar et al.l[200fil : lFedele et al.ll20inl ). We use the model of lGalvet fc Gullbring! \l^m ) 



to predict the amount of FUV emission produced in the accretion column. This model assumes 
that the circumstellar disk is truncated at a few stellar radii by the stellar magnetic field. At 
the truncation radius, gas is channeled along the magnetic fields lines from the disk to the star, 
impacting the stellar surface at approximately the free fall velocity. Such fast moving material 
produces an accretion shock on the surface of the star which heats the gas to ~1 MK. Emission 
from this hot gas heats both the photosphere below the shock and the pre-shock material, producing 
an accretion column spect rum with a significant excess over the photospheric emission in the UV 
(|Galvet fc GullbrinelEggsl ^. 



As 



m 



Siess et al. 



3 uts to the accretion column model, we assume a 0.8Mq source evolving along the 
(j200ol ) evolutionary tracks, from which we obtain the radius, effective temperature, and 
luminosity of the source from 1-10 Myr. We assume the percentage of the stellar surface covered 
by accretion columns, the fillin g factor, remains constan t at 1% and vary the energy flux in the 
accretion column for a given M (jGalvet &: GullbringjI 19981 ) . Using these quantities we calculate the 
column er nission and integra te the flux from 1230-1800 A for the highest and lowest M's observed at 
each age (jGalvet et al.ll2005l ). We also calculated the FUV accr etion emission for the M 's predicted 
from a viscously evolving disk with an initial mass of O.IM© (jHartmann et al.lll998l ). Using this 
technique, we find we can approximate the range of observed FUV luminosities well (left panel of 
Figure [7]), but for typical GTTS masses, the predicted FUV fluxes from the viscously evolving disk 
are too low (right panel of Figure [7]) . 



The missing FUV emission is probably due to several sources which are not included in our 
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accretion column models. First, emission from H2 gas was observed in the FUV spectra of CTTS 
and was necessary to explain the FUV continuum between 1500 and 1700 A. The contribution of 
gas emission can be a significant portion (~ 30%) of the total FUV emission, althou gh it does not 



contr i bute to the disk he a ting, since it likely o r iginates in the upper laye rs of the disk (jlnglebv et al 



2009 : iBergin et all l2004l : iHerczeg et all I2OO2I . l2004l : [France et al.ll2010l ). In addition, we measure 
the flux in the strong FUV atomic lines of existing spectra, which our accretion column models 
do not attempt to reproduce and find that ~ 35% of the total FUV flux comes from hot lines. 
Another contribution necessary for re-producing the short wavelength FUV continuum {^^ 1300 — 
1400 A) could originate in post-shock emission escaping absorption by the pre-shock material. 
The post shock spectrum peaks in the soft X-rays but extends into the FUV around 1300 A 
(jCalvet fc Gullbringill998l ). This emissio n would be import ant for disk photoevaporation models as 
soft X-rays can efficiently heat disk gas (jOwen et al.ll2010l). but may on ly play a role when the M 
drops and the accretion columns become less dense (jGorti et al.ll2009bl ). Additional observations 
are required to determine if this is a viable contribution to the FUV spectrum. 

After accretion has ended and the gas has been dissipated at >10 Myr, FUV luminosities con- 
tinue to decline. This is also observed in magnetically active young stars; all hi gh energy radiatio n 
fields decrease as the sources evolve towards and even along the main sequence (jRibas et al.ll2005l ). 
Given this fact, we looked for any relation between the FUV and X-ray emission in the CTTS and 
non-accretors. In Figure [8] we distinguish non-accretors (diamonds) from CTTS (squares) and indi- 
cate hotter stars by larger symbols. We see no correlation between FUV and X-ray emission in the 
CTTS (with Pearson correlation coefficients of 0.3 and -0.2 for total and normalized luminosities, 
respectively) but we note that the sources with earlier spectral types tend to have the strongest 
total FUV and X-ray emission, as they are observed to occupy the upper right hand corner in the 
left panel of Figure El FUV and X-ray luminosities normalized by the stellar luminosity are indis- 
tinguishable by spectral type (right panel of Figure [8]). We do observe a weak correlation between 
both total and normalized X-ray and FUV luminosities when considering only the non-accreting 
sample, with Pearson correlation coefficients of 0.74 and 0.80, respectively. 

We compare results from our sample of non-accreting young stars (excluding the lo w mass 
TWA sources) to an analysis of solar analogs between 0.1 and 10 Gyr (iRibas et al.ll2005l ). which 
showed that high energy fields between 1 and 1200 A decreased with age, with emission from hotter 
plasmas diminishing more quickly. For our sample of non-accreting solar type sources between 15 
Myr and 1 Gyr, classified as young suns in the Formation and Evolution of Planetary Systems survey 
(FEPS; Meyer et al. 2006), we compare the observed luminosities to tho se predicted a ssumi ng the 
power law relations between fiux in a given energy range and age from I Rib as et al.l (j2005l ). The 
predicted decay of X-ray luminosities is similar to that which we observe for the youngest of the 
solar type non-accretors in our sample but fail to re-produce the 1 Gyr source, HD 53143 (see 
Figure [9]). We also compare the slope of the predicted FUV fiuxes between 920-1180 A to the 
slope of our observations which cover 1230-1800 A. The total predicted luminosities were scaled, 
due to the use of different wavelength bins, but the slope was held constant. Again, the younger 
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sources follow the predicted slope but the 1 Gyr source has an unexpectedly high FUV luminosity. 
It is unclear why HD 53143 is discrepant in both X-rays and FUV f luxes for its age, with several 
different age indicators providing consistent values (jSong et alj 120001 ) . 



Our young solar analogs extend the radiation evolution observed in iRibas et alJ (120051 ) to 
younger ages, and agree with their result that the higher energy r a,diation field s decl ine more 
quickly than the cooler fields. Combining our results with those of iRibas et al.l (|2005l ). we see 
evidence that even past the 10 Myr time period for disk dispersal and planet formation, high 
energy fields strongly in excess of the solar values persist and may have an effect on evolving disks 
and planets. 



4.2. The He II A1640 and C IV A1549 lines 



Using observations for a limited number of T Tauri stars, [Alexander et al.l (|2005l ) have inves- 
tigated the origin of the He II A1640 and C IV A1549 lines and their r elationship to the conti nuum 
flux, responsible for photoevaporation of the gas disk. In particular, [Alexander et al.l (j2005l ) pro- 
pose that the ratio of the He II to C IV is a probe of the EUV radiation field. They base their 
conclusion on the fact that single temperature and density models of collisionally ionized plasmas 
cannot reproduce the observed strengths of the He II line, from which they infer that the line must 
be produced by radiative recombination. Since the ionization threshold for He I is 228 A, they 
conclude that the line is a probe of the strength of the radiation field around this wavelength, 
namely, the EUV. 

Here, we revisit these lines using our larger sample of FUV spectra of T Tauri stars. We 
measure the He II and C IV fluxes in our ACS/SBC and STIS samples of sources with the STIS 
observations degraded to the low resolution of ACS/SBC. As a result of the resolution, we are 
measuring a blend of emission lines in each line feature, including both components of the C IV 
doublet and numerous H2 lines present in the vicinity of C IV and He II. While adding to the 
scatter in our measurements, the dominant line in each region is the He II or the C IV line, with 
H2 line contributions < 10%. As discussed in §3.41 the largest source of uncertainty comes from 
the reddening correction of the FUV spectra, so to estimate errors on the He II and C IV line 
luminosities we assume a range of ±0.2 magnitudes in Ay- Figure [TOl shows He II and C IV line 
luminosities and their dependence on age and Lace- Both He II and C IV line luminosities are 
observed to decrease with age with the lowest values obtained for non-accretors. We find strong 
correlations between line luminosities and Lace, with Pearson coefficients ~ 0.84 and 0.88 for He II 
and C IV , respectively . The stren gth of both lines is clearl y related to accretion, as was previously 
found bv lCalvet etlH and I Johns-Krull et al.l (|20od l. 



In the left panel of Figure [TT] we show the He II to C IV ratio for the CTTS and non-accretors 
of our sample as a function of age. T he value of the ra t io is higher in non-accretors than in the 
CTTS, confirming the trend found by [Alexander et al.l (120051 ). These authors suggest that this 
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trend indicates that the ionizing flux is not powered by accretion, since the ratio would be expected 
to decrease if that were the case. We suggest, instead, that the difference between the He II to C 
IV ratio between the CTTS and the non-accretors arises from the different formation mechanisms 
of the lines in each type of star. The He II A1640 and C IV A1549 lines in non-accretors must 
form in the transition region, as they do in the sun and in magnetically a.ctive stars, indicated 
by th e similarity of the ratio in the two classes of stars (jHartmann et al.l Il979l : iByrne fc Doyle 
19891 ). In contrast, the correlations between line luminosity and Lace found in CTTS indicate that 
these lines must form in the accretion flows or at least in regions powered by accretion energy. 
However, despite the line luminosity dependence on Lacc^ we find no correlation between the He 
II to C IV ratio and Lace, as shown in the right panel of Figure [TTJ Studies of the optical and 
infrared lines of He I and He II in CTTS indicate that the lines have multip le components arising 



in the accretion sho cks, the magnetospheric infall regions, and the winds (jBeristain et al 



2001 



Edwards et al.ll2006l ). No comparable analysis of the C IV line has been done, but it is also likely 
that this line h as contributions fro m the accretion shock and from the infall region, given its large 
velocity width ()Ardila et al.ll2002l ). This diversity may be responsible for the large dispersion of 
the He II to C IV ratios in CTTS. We therefore caution against using the He II to C IV ratio as an 
indicator of EUV emission without a careful treatment of accretion related emission in these lines. 

Regardless of the effectiveness of the He II to C IV line ratio in tracing the E UV emission 



strong EUV emission is likely present in the high energy spectrum of T Tauri stars. iRibas et al 



(j2005l ) showed that EUV emission decreases with age, or increases with the spin of the source. 
They also showed that younger sources have enhanced EUV line emission, in particular strong He 
II and Fe hues around 300 - 350 A. Whether the EUV emission continues to increase at <10 Myr 
like the FUV emission or flattens and remains constant like the X-ray emission is not clear, but in 
either case the EUV emission is expected to be strong. In our model of the accretion column which 
contributes emission to the FUV spectrum ( §4.ip . we find that the EUV luminosity produced in the 



shock is ~ 10^^ ergs~^. In the current model, only ~2% of this emission escapes the high density 
material which is still accreting; however, recent work has shown that varying accretion properties, 
like the accretion r ate or infall velocity, affects the amount of shock emission which can be observed 
(jSacco et al.ll2008l ). Future models of the emission from the accretion column which take this into 
account will provide better estimates of the EUV emission which is available for photoevaporation. 



4.3. High Energy Radiation Fields and Photoevaporation 



In contrast to early photoevaporation models in which the primary mechanism for dispers- 
ing the disk ut ilized EUV radiation f rom the star producing winds with mass loss rates of ~ 
10^^'^ Mq yr~^ ( Hollenbach et al.lll994l ). recent models stress that either FUV or X-ray radiation 
fields are significantly more important in dispersing the disk with mass loss rates 100 times higher 



than in the pure EUV case. While in agreement on the v alues of the mass lo ss rate, I Owen et al 



(j2010l ) achieve these high rates with X-ray radiation while iGorti et al.l (j2009bl ) find the FUV radi- 
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ation drives the photoevaporation. With the large range of ages in our X-ray and FUV samples, 
we explore the photoevaporation theories in the context of our observations. 



Our results agree with previous findings by iPreibisch fc FeigelsonI (j2005l ) and iMercer et al 



(|2009l ) showing that the X-ray luminosity remains roughly c onstant during the first 10 Myr. Re- 
garding models of X-ray photoevaporation. lOwen et alJ ( 20ld ) use Lx = 2x 10^° ergs~^ for hv > 0.1 
keV as their X-ray luminosity. We find similar, though consi s tently lower. X-ray luminosities for 
the typical CTTS mass sources in our sample. lOwen et al.l (|2010l ) note that their input X-ray 
luminosity is large and that further analyses will include lower X-ray luminosities. 

An important new result in this work is the evolution of the FUV radiatio n, where we see that 
the luminosity between 1230 and 1800 A clearly drops between 1 and 10 Myr. iGorti et al.l (l2009bl ) 
incorporate the evolution of FUV luminosity in their models for disk dispersal, assuming that M 
evolves viscously and the FUV luminosity is 0.4 Lace- In Figure [7| we showed the evolution of the 
FUV luminosity along with the predicted FUV emission produced in an accretion column from a 
viscously evolving disk onto a 0.8Mq star. We found that the predicted FUV accretion emission 
is lower than the observed FUV fluxes, although, within an order of magnitude (assuming typical 
mass CTTS have Lace < 0.5 L(7)) and we disc ussed origins for the remaining FUV emission in ^4.11 



The FUV fluxes used by iGorti et al.l (l2009bl ) are similar to those we observe, but we found that 
~l/3 of this emission may be produced by H2 in the circumstellar disk (Ingleby et al. 2011) and 
therefore not directly irradiating the gas in the disk. 

Our observations support the r equir emen ts of both the X-ray a ,nd the FUV dominated photo- 
evaporation models of I Owen et al.l (|2010l ) and lGorti &: HollenbachI (|2009al ): however, several issues 
remain unresolved. For example, circumstellar disks of gas and dust have been observed around 



much older sourc es, including the 10 Myr CTTS T^ 



transitional disk (ICalvet et al.ll2002l : iNa 



around the 7 -10 M yr CTTS in 25 Ori ( 



ita et al. 



lya which maintains a 0.06 Mq disk, albeit a 



Ingleby et al 



201c). Fu l l disks of gas and dust ar e also observed 



2OO9I : iHernandez et al.ll2007al ). According to 



Owen et al.l (j2010l ). the continued existence of full disks at late stages could be explained by the 
range of X-ray luminosities observed, with lower X-ray luminosities increasing the disk lifetime 
and, alternatively, the presence of young WTTS may be due to high X-ray luminosities in those 
sources. In our 25 Ori observations in which ~90% are WTTS, we flnd that the median X-ray 
luminosity (log Lx) is 29.7 erg s^^ and the four accreting sources have only marginally lower or, 
in one case, even higher X-ray luminosities. With similar X-ray luminosities to the WTTS, these 
accreting sources should have cleared their inner gas disks in similar timescales. If we look at the 
X-ray luminosities of the WTTS in the 2.5 Myr Chamaeleon region, we find a median log Lx = 29.0 
s~^, lower than all of the 25 Ori accretors, indicating that if these X-ray luminosities may disperse 
an inner gas disk in '^2.5 Myr, they would clearly be gone by 10 Myr. An additional parameter, 
besides a range of X-ray luminosities, is necessary to explain the diversity of disks observed at a 
given age. 



When including the time evolution of the FUV radiation fields and their dependence on disk 
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properties, like M, iGorti et alj (|2009bl ) find FUV fields evolving on a range of timescales related 
to differences in the accretion properties. Therefore, photoevaporation dominated by evolving 
FUV fields may disperse disks in a range of timescales. With FUV driven photoevaporation and 
evolving FUV fields, the diversity of disks observed at any age may be expected. This appears 
to be a plausible scenario given the large range of observed FUV luminosities and inferred M's. 
Again, the full disks in the 25 Ori association present an obstacle to the theory. The observed FUV 
luminosities of the 25 Ori accretors are ~ 2 orders of mag nitude highe r thaii the WTTS in the 10 
Myr TWA association, and given our results and those of lRibas et al.l (|2005l ). were likely stronger 
at 1 Myr. These FUV luminosities should be high enough to photoevaporate the gas disk, if FUV 
emission was responsible for the gas depletion in the TWA WTTS. Possible explanations for the 
remaining full disks in 25 Ori, in the context of both X-ray and FUV photoevaporation theories, 
i nclude low disk v iscosities or high initial disk masses which increase the time for disk dispersal 
dOwen et al.ll20ld ). 



4.4. FUV Radiation and Dust Evolution 



So far, we have focused on the gas in circumste llar disks; howeyer, it is unclear if gas and dust 



evolut ion proceed concurrently or independently (jPascucci et al.l l2006l ) . [Alexander &: Armitage 
(|2007l ) modeled the affect of photoevaporation on the gas and dust in the circumstellar disk and 
found that the removal of gas from the disk has significant effects on the dust. When photoevapo- 
ration creates a gap in the gas disk, the inner disk dust is also removed. Also, pressure gradients 
formed when gas is removed cause the dust to migrate and may enhance grain growth. For these 
reasons, we may expect to see dust evolution dependent on the strength of the photoevaporative 
fields. With constant X-ray fields throughout disk evolution, we focus on the accretion related FUV 
radiation fields. 



Fedele et al.l (|201Cll ) find that dust dissipation and accretion are characteri zed by different 



times cales and to test this further, we compare FUV luminosities to 2MASS J, H, K (jSkrutskie et al 



20061 ). IRAC {Spitzer InfraRed Array Camera), and MIPS {Spitzer Multiband Imaging Photome- 
ter) colors. IRAC and MIPS photometry, along with disk cla ssifications for o ur source sam- 



ple, w ere taken from the literature: Tau rus (ILuhman et al.l l2010ll. Cham aeleon I (ILuhman et al 



20081'). Chamaeleon II (lAlcala et al.ll2008l^. 25 Ori (IHernandez et al.ll2007al'l. TWA (|Low et al 



2005 



Hartmann et al. 2005 ). and debris disks ( Carpenter et al. 20081 : Chen et al. 2005 ). Figure [12] shows 
the relations between FUV luminosity and inner disk dust tracers. At A < 2 4 fj,m we are tracing 
the dust in the inner regions of the disk; as shown in iD'Alessio et al.l ()2006l ) , 60% of the 24 fim 
emission is coming from inside 10 AU when the dust has not settled towards the midplane. The 
percentage of emission coming from within 10 AU increases when observing at the shorter IRAC 
wavelengths, or when observing a disk in which some degree of dust settling has occurred. 

Among the CTTS (circles) in Figure [121 we find tentative correlations between FUV luminosity 
and each of the disk evolution tracers, with the most convincing correlation in J — K. However, we 
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are likely observing two quantities which correlate with Lace- We showed in Figure[7]that the sources 
with highest FUV luminosities are also the strongest accretors. Disk properties also depend on the 
accretion rate through the disk and the IR spectral energy distributions reflect this; models predict 
that lower M's will produce bluer IR c olors and that the IR spectrum becomes increasingly flatter 
as M increases (jP'Alessio et al.ll2006l : lEspaillatl l2009l ) . Therefore, it is likely that these apparent 
trends are a result of accretion affecting both the FUV luminosities and the tracers of dust in the 
disk instead of a consequence of direct processing of the dust by FUV radiation. This seems to 
indicate that if it is the FUV that is responsible for gas disp ersal, the gas and dust are evolving 
independently. Supporting this conclusion, the transitional (jCalvet et al.l |2002| . l2005l : iKim et al 



20091 ) and pre-transitional (jEspaillat et al.ll2008l . I2OIOI ) disks, the sources likely to have undergone 
dust processing and evolution, have the same FUV luminosities as the full disks, indicating that 
the FUV is not likely having a strong effect on the dust. 



5. Summary and Conclusions 

We present new observations of T Tauri stars in the Chamaeleon I and II and the 25 Ori 
star forming regions. Included are Chandra X-ray and ACS/SBC FUV observations which provide 
constraints on the photoevaporative fluxes present in young stars. Combined with previously 
published X-ray and FUV observations, we explored the evolution of high energy radiation flelds 
with the following results: 

• With the addition of '^40 sources with X-ray observations in the 7-10 Myr 25 Ori star forming 
association, a significant addition to the ~25 previously observed 6-10 Myr T Tauri stars 
(above the substellar limit), we confirm results that there is no decay in the X-ray luminosity 
between 1 and 10 Myr. With consistently strong X-ray radiation fields, disk photoevaporation, 
if driven by X-rays, should proceed on short timescales. However, observations of sources 
with remaining disk gas emission and strong X-ray luminosities at 7-10 Myr and previous 
observations of TW Hya (a source with a transitional disk and strong X-ray luminosity) still 
cannot be explained by X-ray photoevaporation. 

• We show, for the first time, the observed evolution of the FUV radiation fields in pre-main 
sequence stars from 1 Myr to 1 Gyr. We find that the FUV luminosity decline can be 
explained, in part, by the decline of the accretion rate onto the star but note that some FUV 
emission must come from an additional source which is likely gas in the inner disk regions. 
If FUV drives photoevaporation of the gas disk, the diversity of disks observed at 1 Myr 
can be understood by the large range of M's observed at 1 Myr. However, gas disks at 7-10 
Myr require additional explanation, as viscous M evolution predicts that these sources should 
have had strong FUV radiation in the past, so disk dispersal by FUV radiation should have 
proceeded quickly. 

• We show that the power law relations for the decay of high energy radiation fields in solar 
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type active stars described by iRibas et al.l (|2005l ) extend to the significantly younger sample 



discussed here. With both the X-ray luminosities and the slope of the FUV relation accurately 
predicted, we expect that the difficult to observe EUV radiation fields will similarly obey the 
power law trends. This result tells us that strong radiation fields, necessary to remove the 
gas, are present throughout the disk dispersal phase. 

• We find no correlation between FUV luminosities and tracers of dust evolution in the disk 
and furthermore, we find that sources with significant dust processing, the transitional disks, 
have similar FUV luminosities to full disk sources. Therefore, if FUV radiation is the key 
emission field in gas evolution, it is not driving the dust depletion. This would indicate that 
gas and dust evolution occur independently, though interestingly on similar timescales. 
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-5 






0.066 




11, 12 


OBla 1626 




K8 


4.1 


X 


10- 


-4 






>0.83 




11, 12 


OBla 1630 




M2 


6.6 


X 


10- 


-5 


1.0 X 10-3 


0.16 


0.09 


11, 12 


OBla 1663 




M5 


4.1 


X 


10- 


-5 






0.10 




11, 12 


OBla 1755 




M3 


1.0 


X 


10- 


-4 






0.21 




11, 12 
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Table 1 — Continued 



Source Alt. Name 


SpT 


Lx 




Lpuv 






Lace 


Reference"''' 






















RX J0526.7+0143 


KO 


1.0 X 10" 


-3 


- 




- 


- 




1 A 
14 


J05244498+0159465 


K6 


8.3 X 10" 


-4 


- 




0.66 


- 




1 1 
il 


J05260639+0137116 


M2 


6.6 X 10" 


-5 


- 




0.041 


- 




1 1 


J05263125+0141492 


M4 


5.2 X 10" 


-5 


- 




0.033 


- 




1 1 


TWA 




















TW Hya 


K7 


6.6 X 10- 


-4 


1.3 X 10" 


-3 


0.32 


0.03 


1 


lb 


HD 98800 


K5 


1.6 X 10' 


-3 


2.6 X 10" 


-5 


0.66 





1 


16 


TWA 7 


Ml 


1.0 X 10- 


-4 


1.0 X 10" 


-5 


0.32 





1 


16 


TWA 13 A 


Ml 


2.6 X IQ- 


4 


1.6 X 10" 


o 


0.16 





1 


16 


TWA 13B 


Ml 


2.6 X IQ- 


-4 


1.6 X 10" 


-5 


0.16 





1 


16 


Field 




















HD 12039 


G4 


1.6 X 10' 


-4 


4.1 X 10" 


-5 


0.33 





1 


16 


HD 202917 


G5 


5.2 X 10- 


-4 


6.6 X 10" 


-5 


0.66 





1 


16 


HD 53143 


Kl 


2.6 X 10- 


-5 


8.3 X 10" 


-6 


0.66 





7 


16 


HD 61005 


G8 


4.1 X 10- 


-5 


1.6 X 10" 


-5 


0.66 





1 


19 


HD 92945 


K2 


2.1 X 10- 


-5 


8.3 X 10" 


-6 


0.41 





1 


16 


MML 28 


K2 


5.2 X 10- 


-4 


3.3 X 10" 


-5 


0.41 





1 


16 


MML 36 


K5 


8.3 X 10- 


-4 


1.0 X 10" 


-4 


1.0 





1 


16 



Note. — 

Values in italics were taken from the literature while those in regular font were derived in this work. 
Sources in bold were observed with STIS rather than ACS/SBC.__ 
"Spectral t ype, stell a r, and acc retion luminosit y re ferences. 1 



therein, 2) 



Luhman 



Briceho et al 



Espaillat et al. tOOS). 7) iKalas et al.l (|200i 



(120041). 3) iLuhmanI (l2007l l. 4) ISpezzi et al 



Ingleby et a. 



Nguven et al 



(|l994l ). 14) lAlcala et al.l hmd ) 



(120081 ). 5) 



200! 



(120071). 11 ) Brice fio et a. (2010), 12) iHernandez et al.l (l2007al ). 13) iHuenemoerder et al 



IZ9) 



(120091) and refere nces 



Calvet et a. 



Briceho et al 



(I2OO5I ). 6) 
(I2OO5I). 10) 



''X- rav luminosity referen ces. 1 5)lGiidel et al.l tom. 1 6) IVoges et al.l (119991 ) . 17) IVoges et al.l (j200nl ) 
18) IPeigelson et aD (|l993l ). 19) Iwichmann et al.l (|2003l ) 



LaccS for T26 and T54 were calculated from the U band excess, whereas the remaining Chamaeleon 
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sources have Lace calculated using Ha 10% equivalent widths observed here or taken from the literature. 
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Table 2. X-ray Detections: Possible and Confirmed T Tauri Members 



2MASS Source 


SpT 


EW Ha 


EW Li 6706 A 


log Lx 


05243490+0154207 


K4-K5 


-0.46 


0.10 


29.0 


05244265+0154116 


G9-K0 


1.54 


0.32 




05244320+0200355 










05251394+0143313 










05265823+0136078 










05270173+0139157 


K6 


-0.39 


0.17 


29.3 


10575375-7724495 











Note. — Spectral types were determined using the Ca I absorption 
lines. With the early spectral type, additional membership confir- 
mation is needed for 2MASS 05244265+0154116 as Li is not a good 
indicator. 



Table 3. X-ray Spectral Analysis 



Source 


Net Counts 


Nh 




kTa 


^red 


Flux 








10^1 cm-2 


keV 


keV 




erg s ^ cm ^ 


CR Cha 


1996 


0.8 


0.84 


2.97 


0.95 


1.4x10" 


-12 


CHXR 30A 


324 


12.2 




2.14 


1.01 


4.6x10" 


-13 


SY Cha 


283 


0.8 


0.34 


2.94 


0.81 


2.4x10" 


-13 


SZ Cha 


692 


0.8 


0.85 


4.53 


0.74 


3.7x10" 


-13 


TW Cha 


313 


0.8 


0.41 


1.95 


1.08 


1.8x10" 


-13 


CVSO 38 


266 


1.1 


0.98 


6.98 


1.42 


4.1x10" 


-13 


J05244498+0159465 


285 


1.1 


0.32 


1.59 


2.06 


2.5x10" 


-13 


RX J0526+0143 


369 


1.1 


0.31 


1.48 


2.29 


2.6x10" 


-13 
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Fig. 1. — H-R diagram of X-ray and FUV samples of Chamaeleon sources. We show the locations 
on the H-R diagram of Chamaeleon I X-ray and FUV sources and identify those with both X-ray 
and FUV observations. All WTTS were observed in X-rays and do not have FUV spectra. We 
also show the location of the older Ch amaeleon II sourc es which were observed with ACS/SBC. 
Solid lines are evolutionary tracks from lsiess et al.l (200C) and are labeled in units of solar masses. 
Dashed lines are isochrones and correspond to log age = 5.5, 6, 6.5, 7, and 7.5 years. 
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Fig. 2. — H-R diagram of 25 Ori association X-ray detections. We show the location of all 25 Ori 
X-ray sources on the H-R diagram with evolutionary tracks and isochrones as defined in Figured) 
Sources fall primarily between the 6.5 and 7 My r isochrones with earlier (hotter) sources appearing 
older due to birth line effects (iHartmannl l2003l ) . 
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Fig. 3. — Ha line profiles from MIKE. We show the accreting sources in our sample which have 
MIKE observatio ns, ordered b y log M (listed below each source name) calculated from the 10% line 
width relation of lNatta et all koo4 ). The 7-10 Myr sources, CVSO 206, OBla 1630, and CVSO 
35, have remarkably wide and asymmet ric Ha profiles, ind icating that accretion is strong and gas 
remains in the inner circumstellar disk ([inglebv et al.ll2009l ). 
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Fig. 4. — Chandra ACIS spectrum of the accretor CR Cha. We show the extracted spectrum of 
CR Cha (the source with the highest counts in our sample) along with our best APEC model fit 
to the spectrum. 
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- — 25 Ori 

- - - Cha I 

6 - 



z 4 - 

- r 

2 
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28.5 29.0 29.5 30.0 30.5 

log X-ray Luminosity (erg/s) 



Fig. 5. — Distributions of X-ray luminosities for 25 Ori and Chamaeleon I samples. The distribution 
of X-ray luminosities, corrected for absorption, in each sample does not reproduce th e lognormal 
shape of the X-ray luminosity function found for COUP sources in lFeigelson et al.l (j2005l ). This indi- 
cates that our samples are likely incomplete in mass. For this reason we supplement our Chamaeleon 
I observations with previous X-ray observations; however, no additional X-ray observations exist 
for 25 Ori. 
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Fig. 6. — Left: Median X-ray luminosities of young star forming associations vs age. For each 
association, we plot the median and the full range in X-ray luminosity for sources with spectral 
types K5 and later. The diamonds represent X-ray lun iinosities calculat ed in thi s work and the 
boxe s represent thos e taken from the lite rature: T aurus (Giidel et al. 20071'). C ha I ( Feigelson et al 



199S 



1993) . Tru mpler 37 (iMercer et al.l 12009). r? Cha (ILoDez-Santiago et al.ll20ld ). TWA (|Voges et al 



2000l ). We include the sample of iFeigelson et al.l (|l993l ) for Cha I (offset in age for clarity) 



because with only two fields observed here, the Cha I sample was limited. Our Cha I observations 
reproduce the range observed by ROSAT well and actually probe dimmer X-ray sources. Right: 
Normalized X-ray luminosity vs age. Here we divide the observed X-ray luminosities by the stellar 
luminosities to normalize by spectral type. Both figures show that between 1 and 10 Myr, there is 
no decline in X-ray luminosity. 
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Fig. 7. — Left: FUV luminosity vs age for young low-mass stars from the T Tauri phase up to 1 
Gyr. The FUV luminosity of CTTS, displayed as asterisks when accretion rates are known and 
diamonds when they are not, is observed to decrease with age until becoming entirely chromospheric 
in WTTS (triangles). Upper limits were calculated for non-detections, primarily sources with high 
Ay, and are represented by downward arrows. Solid lii ies show the predicted F UV luminosity 
produced in an accretion shock onto the stellar surface (jCalvet Gullbring||l998l ) as the source 
evolves and M decreases. The top and bottom lines are calculated for the observed range in M 
at each age and the middle line is calculated fr om M predictions in which a 0.8 Mp) CT TS with 
an initial disk mass of 0.1 Mq evolves viscously (jHartmann et al.lll998l : ICalvet et al.ll2005l ). Right: 
FUV luminosity vs Lace- The solid line shows the predictions of FUV luminosity from accretion 
onto the 0.8 Mq CTTS through the viscously evolving disk. Here, we omit the predicted FUV 
luminosities for the upper and lower limits of observed M's, which were shown on the left. 
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Log X-roy Luminosity (erg/s) Log X-roy Luminosity / Stellor Luminosity 

Fig. 8. — Left: Comparison of total FUV and X-ray luminosities. Boxes represent accreting sources 
and diamonds represent non-accreting sources. The largest symbols are G stars, medium symbols 
are K stars and the smallest symbols are M stars. There is no correlation between FUV and X-ray 
luminosities among this sample; however, the non-accretors have the lowest FUV luminosities. The 
subset of non-accretors has weakly correlated total X-ray and FUV luminosities. Right: Comparison 
of fractional X-ray and FUV luminosities. Again, no correlation is observed but non-accretors have 
the lowest fractional FUV luminosities. One CTTS, CR Cha, appears in the region populated by 
non-accretors (< —3.5). We note that this is the only Chamaeleon I source with Ay = 0. 
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Fig. 9. — The luminosity of young; solar ana logs compared to power law relations of high energy 
radiation fields derived in iRibas et al.l ( 2005 ). The solid line s hows the luminosit ies predicted for 
the summed 1- 20 A and the 20-100 A flux bins as defined bv lRibas et al.l ( 20051 ). approximately 
the same wavelength region as our 0.2-10 keV observations. The dashed line shows the slope of the 
predicted fluxes from the 920-1180 A wavelength bin but scaled vertically in luminosity to fit the 
FUV luminosities of our longer wavelength observations. The slope of the FUV power law relation 
provides a good fit to our observed FUV luminosities. The predictions fail to fit the 1 Gyr source 
in both the X-ray and FUV regimes. 
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Fig. 10.— Top Left: He II A1640 A luminosity vs age. Bottom Left: C IV A1549 A line luminosity 
vs age. Symbols are defined as in Figure [71 In both cases the luminosity of the line is observed to 
decrease with age. Top Right: He II A1640 A luminosity vs Lace- Bottom Right: C IV A1549 A 
line luminosity vs Lace- Both He II and C IV are correlated with Lace- 
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Fig. 11. — Ratio of He II to C IV as a function of age (left) and accretion luminosity (right). 
Symbols are defined as in Figure [71 We see that the non-accretors clearly have a higher value of He 
II/ C IV, approximately 1, while the CTTS have He II/ C IV < 1; however, we see no clear trend 
with accretion luminosity. 
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Fig. 12. — FUV radiation and tracers of disk evolution. We compare the FUV luminosities to several 
indicators of evolution in the disk. [3.6], [8], and [24] represent IRAC and MIPS photometry at 
3.6, 8, and 24 /xm respectively. Open symbols are sources with spectral types earlier than K5 
and closed symbols are spectral types K5 and later. J-[3.6], A'-[8], and i^-[24] probe regions 
progressively further in the disk. We see only tentative correlations between high energy emission 
and dust evolution tracers, but this is likely due to both quantities being dependent on accretion 
rate. We do see that the full and transitional or pre-transitional disks (boxed symbols) have similar 
FUV luminosities, indicating that the FUV radiation is not very significant in dust evolution. 



